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Abstract 

We resum a class of large Sudakov logarithms affecting Higgs boson production in the exclusive 
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next-to-leading order cross section and present a detailed numerical study assuming realistic LHC 
cuts. Careful attention is paid to the matching procedure, and to the theoretical uncertainties 
induced by residual scale variation. We find that the matched NLL' -|- NLO cross section has 
significantly smaller uncertainties than the fixed-order result, and can be used to alleviate the 
theoretical errors hindering current Higgs analyses at the LHC. 
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I. INTRODUCTION 



The discovery last year of a new boson by the ATLAS and CMS collaborations at the 
LHC [H 12] has ushered in a new era in particle physics. The future program of the LHC, 
and the next stage of experimental studies in high energy physics, will be largely devoted 
to measuring and understanding the properties of the new state in order to determine the 
underlying theory from which it arises. The initial data provides only a hazy glimpse at the 
properties of the new particle. Initial measurements of its branching ratios into various final 
states indicates that its couplings are consistent with those predicted for the Standard- Mo del 
Higgs boson [^ij, as is its parity Significant work will clearly be needed to sharpen our 
picture of the new state. 

A major component of the quest for a better understanding of the newly-discovered state 
will be the improvement of theoretical predictions for Higgs-boson production and decay 
channels. It is well-known that the next-to-leading order (NLO) QCD radiative corrections 
to the gluon-initiated gg ~> H + X production process are so large [5HTT] that next-to-next- 
to-leading order (NNLO) results are required to realistically describe LHC measurements. 
NNLO calculations are available for inclusive Higgs boson production [T2HI1] and for the 
fully-differential Higgs plus zero-jet cross section [T5HI9]. Recently, first results for the NNLO 
cross section in the Higgs plus one-jet channel have become available [2Dj. As the gluon- 
fusion process is the primary production mode at the LHC, it has received an enormous 
amount of theoretical attention. In addition to the QCD corrections through NNLO, the 
leading electroweak [211 122] and mixed QCD-electroweak [23] corrections are known. More 
complete recent reviews of precision predictions for the gluon-fusion and other production 
channels can be found in Refs. [^H - I^Tj . 

Most of the theoretical predictions described above are obtained using fixed-order pertur- 
bation theory, and assume that there are no severe cuts on the phase space of the hadronic 
radiation produced in association with the Higgs. Unfortunately, such constraints are present 
in several important Higgs search channels. A well-known example is that of a Higgs boson 
decaying to l^-bosons [28] [29] . The background composition to this signal chang es as a 
function of jet multiplicity. In the zero-jet bin the background is dominated by continuum 
WW production, while in the one-jet and two-jet bins, top-pair production becomes increas- 
ingly important. The optimization of this search requires cuts dependent on the number 
of jets observed, and therefore also on theoretical predictions for exclusive jet multiplici- 
ties. Theoretical predictions for exclusive jet bins suffer from large logarithms of the form 
L = ln{Q / p1^^°) , where Q ~ Mh denotes the hard scale in the process. For disparate scales 
Q and p1p^°, these logarithms can overcome the suppression that occurs at each order 
in perturbation theory, and fixed-order results can consequently lead to incorrect conclu- 
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sions. For example, for the experimentally relevant values p^^° ~ 25 — 30 GeV, residual 
scale variations in fixed-order calculations lead to estimated errors that do not accurately re- 
flect uncalculated higher-order corrections [301432] . The importance of controlling these large 
logarithms in order to obtain reliable central values and uncertainties for the gluon-fusion 
channel in exclusive jet bins has been emphasized in the literature [3214M] . 

The theoretical community has invested significant recent effort in attempting to resum 
jet-veto logarithms to all orders in perturbation theory in order to more accurately model the 
LHC Higgs signal. The resummation for the zero-jet bin cross section in the presence of the 
anti-/cT algorithm was first obtained at next-to-leading logarithmic (NLL) accuracy [30] and 
later extended to NNLL accuracy [35], [36] using two different theoretical approaches. The 
importance of potentially large Ini? corrections on numerical predictions, where R is the jet- 
radius parameter in the anti-/cr algorithm, was studied in Ref. [37] . A significant reduction 
of the residual theoretical uncertainties was obtained in the zero-jet bin by resumming the 
jet-veto logarithms. Given that the theoretical uncertainties are currently one of the largest 
systematic errors affecting the one-jet bin analyses of the Higgs-like particle properties [38] . 
it is desirable to formulate the resummation when final-state jets are also present. 

In a previous paper we established the formalism necessary for resummation of the Higgs 
plus jet process by deriving a factorization theorem using soft-collinear effective theory 
(SCET) [391 - H3] for the production of color-neutral particle and one or more jets in the 
presence of a jet- veto [S]. This result assumes that the transverse momenta of all hard jets 
are larger than the veto scale. We calculated contributions through next-to-leading order in 
the exponent of the Sudakov form factor, and presented initial numerical results for Higgs 
production in association with a jet at the LHC. We found that resummation of the jet-veto 
logarithms significantly improves the reliability of the perturbative expansion, and could 
potentially lead to a reduced theoretical systematic error in experimental studies. 

In this manuscript we extend the calculation of Ref. [Hj in several ways. We first present 
a calculation of the NLO soft function appearing in the factorization theorem for exclusive 
Higgs plus one-jet production. This allows the extension of the resummation accuracy to 
NLL' accuracy, using the logarithmic counting established in Ref. [33]. This level of loga- 
rithmic accuracy implies that we correctly obtain the first three logarithmic corrections at 
each order in the QCD coupling constant: a^L^, UgL and a^; a^L"^, a^L^, and a^L^; a^L^, 
a^L^, ct^L^] and so on. We match our results to fixed-order to obtain a NLL' -|- NLO pre- 
diction, and present numerical results for use in LHC analyses. We first demonstrate that 
the region of phase space where the leading-jet transverse momentum is of order the Higgs 
mass accounts for nearly half of the error in the fixed-order NLO prediction for Higgs plus 
one jet, and is therefore a prime candidate for an improved theoretical treatment. We then 
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perform a detailed study of the residual theoretical uncertainties using our resummed pre- 
diction that accounts for the variation of all unphysical scales remaining in the prediction. 
Even with a very conservative treatment of the errors, a significant reduction of the residual 
uncertainty as compared to the fixed-order estimate is found; the estimated uncertainties 
decrease by up to a quarter of their initial values. Our results, and the improvements in the 
zero-jet bin obtained previously, should form the basis for future theoretical error estimates 
in experimental analyses of Higgs properties. 

Our paper is organized as follows. We review the factorization theorem of Ref. in 



Section un We discuss the extension of the resummation to the NLL' level in Section III 



and present the calculation of the previously unknown one-loop soft function. A detailed 
discussion of numerical results for the LHC is given in Section IV We describe there how 
we estimate theoretical uncertainties in both the fixed-order and resummed results, and 
demonstrate that the resummation of jet-veto logarithms reduces the theoretical systematic 
error affecting LHC analyses. Finally, we conclude in Section |Vj Many technical details 
needed for the numerical studies are given in the Appendix. 



II. REVIEW OF THE FACTORIZATION THEOREM 



We begin by reviewing the salient features of the factorization theorem for exclusive 
Higgs plus one-jet production |31]. The factorization of the cross section into separate hard, 
soft, and coUinear sectors is complicated by the presence of the jet algorithm needed to 
obtain an infrared-safe observable. Following the experimental analyses, we use the anti-/cT 
algorithm |l5] to define jets. Anti-fc-r jets are built using the following distance metrics: 

Pij = mm{p:^^^,p:^j)ARij/R, 

Pi = VtI- (1) 

The anti-/cT algorithm merges particles i and j to form a new particle by adding their four- 
momenta if Pij is the smallest among all the metrics. Otherwise, i or j is promoted to 
a jet depending on whether pi or pj is smaller, and removed from the set of considered 
particles. This procedure is repeated until all particles are grouped into jets. We note that 
AR^j = Aiiij + ^(fiji where Ar^jj and A0jj are the rapidity and azimuthal angle difference 
between particles i and j, respectively. R is the jet-radius parameter, which in practice is 
chosen to be around 0.4 — 0.5. 

We demand that the final state contain only a sing le jet with > ~ 25 - 30 

GeV. Other jets with a transverse momentum above this threshold are vetoed. Since p^*° is 
usually substantially lower than the partonic center-of-mass energy (A = p1^^° / y/^ <ti 1), the 
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vetoed observables are usually very sensitive to soft and coUinear emissions. We will make 
the following assumptions in order to proceed in our analysis: 

pJ^^mn^^s] l>i?2>A2; ^log2/?<l. (2) 

The first assumption leaves us with a two-scale problem and allows the measured final- 
state jet to be described by a separate coUinear sector. The second of these requirements is 
necessary to insure that the measurement function factorizes into separate measurements in 
each of the coUinear sectors. The third requirement ensures that logarithms associated with 
the anti-ZcT parameter R need not be resummed. We will see later that the first assumption 
is satisfied in approximately 30% of the relevant phase space for Higgs plus jet production 
at the LHC, and that this parameter region contributes roughly half of the total error. We 
will therefore be able to improve the theoretical description of a significant fraction of the 
LHC Higgs signal. Given that pV^^° ^ 25 — 30 GeV and R ~ 0.4 — 0.5 , when the leading jet 
~ niH, the second two assumptions are also justified. 

Our effective theory consists of the following low-energy degrees of freedom: 

• a coUinear jet mode with momentum pj = + kj, where nj is the light-cone vector 
along the jet direction; 

• two coUinear modes propagating along the beam axes a and b, with pi = + ki for 

i = a,b; 

• a soft mode with momentum ks- 

The residual momenta fcj, fcj and the soft momentum kg all scale as a/SA, while the large 
components of the three coUinear momenta scale as uJi ~ \/^. Momenta with smaller scalings, 
such as ultrasoft modes, do not contribute to the final-state observable and can be integrated 
over, and therefore need not be introduced. We are able to utilize an effective-theory frame- 
work because of how the anti-ZcT algorithm clusters the soft and coUinear modes. Referring 
to the metrics defined in Eq. ([T]), we find 

pjj ^ Pj pjs ~ R'^ , pja ~ pjb ~ R~^ log , 

Pss ~ Paa ~ Pbb ~ (Ai?)"^ , Psa ~ Psb ~ Pab ~ (Ai?)"^ log A~^ , 

Ps ^ Pa ^ Pb ^ A~^ . (3) 

From pjj and pj in the first line, we see that the initial clustering combines the final- 
state hard emissions into a jet, so that the soft radiation sees only the jet direction and 
does not probe its internal structure. Also, since the clustering between the soft and jet 
radiation typically occurs earlier than the clustering among the soft radiation, clustering of 
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soft particles across the jet boundary is unlikely to happen |l5l HE]. We see from the second 
line that the mixing between the soft and beam sectors is power-suppressed, as is the mixing 
between the beam and jet sectors. Denoting the measurement function that imposes the jet 
clustering and vetoing as jCi, these factors imply that we can factor jCi into the product of 
measurement functions acting separately on the soft, jet, and beam sectors, 

M = MjMsMaMb, (4) 

up to power-suppressed corrections in and R. 

The remaining steps in the derivation of the factorization theorem utilize the standard 
SCET machinery, and are presented in detail in Ref. [H]. The final result for the cross 
section for exclusive Higgs plus one-jet production takes the following form: 

dcTNLL' = d$//d$j J^i<^H, $j) j da^adxfe^ {2Tif5^ {qa + qh-q.j- qn) 

a,b 

^EE^K^ • S)Xa,.,. ® fM^b,,j, ® hix,)Jj{R) . (5) 

spin color 

We have denoted explicitly by the subscript that we will evaluate this cross section to the 
NLL' level, d^n and d^j- are the phase space measures for the Higgs and the massless jet 
J, respectively. $j) includes all additional phase-space cuts other than the pt veto 

acting on the Higgs boson and the hard jet. H is the hard function that comes from matching 
full QCD onto SCET, and S describes soft final-state emissions. The trace is over the color 
indices. The functions X and J describe collinear emissions along the beam axes and along 
the final-state jet direction, respectively. The measured jet should be much larger than 
p^^°. Operator definitions for all functions are given in Ref. As our purpose here is 

to only briefly review the factorization theorem before presenting new results, we do not 
reproduce these definitions explicitly, and instead refer the reader to the quoted reference. 
The tree-level and one-loop expressions for the jet and beam functions needed for numerical 
studies are presented in the Appendix. 

III. EXTENSION TO NLL' 

A primary goal of this manuscript is to extend the resummation of jet- veto logarithms to 
the NLL' level, following the notation of Ref. [33] . This level of logarithmic accuracy implies 
that we correctly obtain the following towers of logarithms: a<jL^, a^L and a^; a^L^, a^L^, 
and a^L"^; alL^, a^L^, a^L^; and so on. The following ingredients are required to obtain 
this accuracy: 
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• the two-loop cusp and one-loop non-cusp anomalous dimensions which control the 
evolution of the beam, jet, soft and hard functions; 

• the one-loop hard, beam and jet functions; 

• the one-loop soft function. 

The requisite anomalous dimensions, as well as the one-loop jet and beam functions, were 
obtained in Ref. They are included in the Appendix for completeness, as is a detailed 
discussion of their implementation into Eq. ([s]). The one-loop hard function for the gg, qg 
and qq channels can be obtained from the literature |17] . The previously unknown quantity 
is the one-loop soft function. Its calculation is rendered non-trivial by the presence of the 
final-state jet. We describe our computation of the soft function below. 



A. Calculation of the one-loop soft function 

We begin by defining the measurement function for the soft sector: 



where we have set Qa,b = 0{a — b) and ARkj = ^/Ay^~+~K(p. The first term allows soft 
emissions that are well-separated from the final-state jet but have a transverse momentum 
softer than p"!^*^", while the second term allows harder emissions that are within the final-state 
jet radius. As the soft function is simply the square of the soft current integrated over the 
allowed phase space, we can immediately write the NLO contribution to the soft function as 

S = -yMzrVT^-T, [ d'kS{k') Ms, (7) 



where the sum is over the two beam directions and the final-state jet direction. The rij 
denote light-like vectors in each of these directions, while the Tj denote color operators in 
either the fundamental or adjoint representations, depending on whether i denotes a quark 
or a gluon. k is the momentum of the gluon emitted from the eikonal lines. We note that 
the leading-order soft function is normalized to unity. Parameterizing 

k^ = /ct (cosh cos 0, sin 0, sinhy) , 

n'^j = (cosh yj, 1,0, sinhyj), (8) 

and setting 

^ARkj,R = 1 - ^RARkj ) (9) 
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we find 



i<j ^ 



2(27r)'^-i ^ ' 'J ^ ' ' ' m-krij-k 

X (^Op^eto + BiJ^AiJfcjQfcr.p^s*" ) (10) 

where fid denotes the d-dimensional sohd angle and = sin (j). 

We will proceed by reducing these integrals as far as possible analytically, although we 
will end up with two remaining integrals which we evaluate numerically. Since the first such 
integral is a pure number, while the second depends only upon R, this does not affect the 
speed of the numerical program we construct. We note that the soft function will have 



rapidity divergences arising from the rii ■ k in the denominator of Eq. (10), which we will 
regulate by multiplying the integrand by the following factor [IS] : 

|2A:,,3|-V''. (11) 

This regulates the rapidity divergence as a pole in rj, which is then removed by renormaliza- 



tion. There are two distinct structures to consider in Eq. (10): the first when both i and j 
denote a beam direction, and the second when i denotes the jet direction. We write the soft 
function as 

S = Ta ■ Tf, Snn + Ta ■ Tj SnJ + Tf, ■ Tj Sfij (12) 

and study each structure separately. We note before continuing that the virtual correc- 
tions are scaleless, and have the effect of converting the infrared poles in the real-emission 
corrections into ultraviolet poles. 

Snn'- This case will have rapidity divergences as n ■ A; — )■ and n ■ — > 0. We can therefore 



replace the regulator of Eq. (11) using 



|2A;,,3|-V^ k-^u^eM-V \y\) ■ (13) 

This relation is valid in the large y limit in which the rapidity divergence occurs. Since we 
are treating 0{R) contributions as power-suppressed terms and there are no singularities 
associated with emissions along the final-state jet direction, we can expand QARkj,R = 1 + 
0{R) to derive 

S^-^'^^C^'^^'i^l (14) 



(27r)"'-i 2e + ?7 V jj, J \ii J r] 

This can be easily expanded in both rj and e in order to isolate the poles, which are then 
removed by renormalization. 
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Snnj'- The calculation of this structure is more involved than the previous one. We first note 
that the n ■ A; — )■ rapidity divergence is correctly regulated by the following replacement, 
valid in the large-y limit in which the divergence occurs: 

|2fcg,3|-V = 2-^-^1 sinhyl-" k-''u''exp{-r]ye{y)) . (15) 



Performing the kq- integration in Eq. (10), and afterwards using the relation Qar^:j,r 
1 - ©R.A/jfej, we find 

_ 2 ^i-2e vr |^z/exp(-yj)^ 

i^nnj — 9s~, 



{2ttY-^ 2e + r] \ ij, J \ /i 
X / dAy / —{s^) — Qar,j,r, 16 

J-oo Jo (cosh Ay - COS 0) 

where we have shifted the integration variables to Ay = y — yj. We can no longer directly 
expand the integrand for i? — )■ 0, as there is a logarithmic divergence in this limit associated 
with collinear emissions along the final-state jet direction. We will proceed by isolating the 
small- i? behavior of the integrand by defining 



X = (S0) 



-2e exp(A|/) exp{-r]Ay QjAy)) 

(cosh Ay — cos (p) 
_2,2exp(-r/Aye(Ay)) 



We divide the integrand into the two structure X — and Xr. The first piece is finite as 
i? — )■ and can be Taylor-expanded in that limit, while the second contains Ini? behavior 
and cannot. We further find it convenient for the X — Xr term to divide the Ay integral into 
two regions: Ay < 0, for which the rapidity divergence cannot occur and consequently we 
can set f] = 0; Ay > 0, in which the rapidity divergence can occur and r] must be retained. 
Performing the indicated manipulations, we reduce Snnj into the sum of three structures: 



ri qlnR I qrj , qri=0 



where we have identified 



Snn. - 9s Ye [-y) L Jo ~ ^'^'^ ' 

q, ^ 2^2-2e 1 fp^'Y''-" fuexp{-yj)y2 

'~'nn T 9, 



(27,)d-i 2e + r] \ fi J \ /x J r] ' 

2e 



^ i^^^ (^^y ' f/-f [X-X«-2(..,)-e(A,)] .(19) 
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We will numerically calculate the integrals which are left upon expansion of the integrands in 
e. We note that the Snnj term of Eq. (12) can be obtained by taking yj — )• —yj in Eq. (19), 
which affects only the S^^^ structure above. 

Combining all of the information presented above, expanding in both r] and e and removing 
the poles via renormalization, we are left with the following final result for the soft function: 



S 



as_ 



veto 

i:2 + 41n^L 



2rjLlni?2 + AyjL{T, ■ Tj - T, ■ Tj 



- iT^ + m^+Tj[c + fiR)] 
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where we have abbreviated L 
by the following integrals: 

c = 4 



(20) 



ln(/i^/p^'^*"'^). The constant c and function f{R) are given 



TT 



dAy 

) Jo 

-4 log(2) log R^ +8 



dAy 



" d0 log(s^) 



TT AR 



e 



kj 



(21) 



As stated above, we determine these quantities numerically in our analysis. As we only 
perform our analysis for a very small set of R values, this is sufficient for our purposes. We 
note that the and L terms agree with those predicted in Ref. [Hj using the cancellation 
of the combined running of the hard, jet, beam and soft functions. For the Higgs production 
process considered in this work, we have the following color identities: for the ggg channel, 
Ti-Tj = -Ca/^; for the gig2fi'3 channel, Ti-T2 = -(Cf-Ca/2) and Ti-Ts = T2-Ts = -Ca/^. 



IV. NUMERICAL RESULTS 



We now present numerical results for exclusive Higgs plus one-jet production at the LHC. 
We first discuss the matching of the resummed result with the NLO cross section to obtain 
a renormalization-group (RG) improved NLL' + NLO prediction, and demonstrate that we 
correctly capture the large logarithms associated with p^^°. We also discuss the parameter 
region in which our effective-theory framework is valid. Although this region, with p'^ ~ rriH, 
makes up only ~25-30% of the signal, it accounts for nearly half of the theoretical uncertainty 
in the one-jet bin. We then describe in detail how we estimate the theoretical uncertainties in 
both the fixed-order and RG-improved results. For the fixed-order cross section we follow the 
"ST" recommendations of Ref. [32] . Our treatment of the theoretical uncertainties of the RG- 
improved result is necessarily more involved, as our effective-theory approach improves the 
prediction over only part of the relevant phase space. We adopt a combination of direct scale 
variation, which is standard in resummed calculation p2] , and the ST recommendations for 
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the fixed-order region, as described below. Finally, we show that the resummation of the jet- 
veto logarithms leads to a sizable reduction of the exclusive Higgs plus one-jet uncertainties 
at the LHC. 



A. Matching NLL' with NLO 

We begin by matching our resummed expression with the fixed-order NLO result to obtain 
a NLL' + NLO prediction. We use the NLO predictions for Higgs plus one-jet contained in 
MCFM [50]. We obtain our prediction by setting 

In this equation, ctnlo is the fixed-order NLO cross section obtained from MCFM, and ct^ll 
is the resummed cross section up to NLL' accuracy presented in Eq. ([s]). ct^'^l' captures 
the singular features of ctnlo? and is obtained by expanding ctnll' with all scales set to a 
common value /x. Schematically, we have 



o"NLL' = ctlo (1 + «.5'o) e ^^^^ 



(asi)-9NLL(asi) 



^NLL' = ^LO (l + as [-92L'^ - 9iL +9o]) , (23) 

where L^'ll and gjsiLL resum the leading and next-to-leading logarithms, respectively. The 
difference between ctnlo and the expanded NLL' result c^^i^, only contains power-suppressed 
contributions for large values of Q: 

(p-veto p!^t° \ 
^-L , logi?, ■ ■ ■ J , (24) 

with L = log{Q/p^^°), and Q stands for any kinematic quantity of order mn- Since the 
scale QR is used to define the jet mode, the R^L terms are regarded as power suppressed. 

To demonstrate that our formalism correctly captures the singular terms at NLO as 
L — )■ 0, we plot in Fig. [ijthe fractional difference between the expanded cross section c^'^l' 
and the NLO result as a function of pl^*^". We note that due to the power-suppressed R^L 
terms, the difference does not completely vanish as p^^° goes to for fixed R. However, we 
see from the several R values plotted in Fig. [T] that in the i? — > limit, these terms also 
vanish, as expected. 

B. Validity of the effective theory 

Having demonstrated that our NLL' result correctly captures the logarithms of p^^*", we 
comment here briefly on its expected range of validity. In our derivation of the factorization 
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FIG. 1: Shown are the fractional differences ((Tnlo — cr^j^j^, ) /cjnlo for R = 0.4 (solid blue), 
0.1 (dashed orange) and 0.06 (dotted red), respectively. The colored bands represent the esti- 
mated numerical uncertainties. The differences between the expanded NLL' and the fixed-order 
NLO calculations are small compared with the total cross section. 



a. 



< 




FIG. 2: Aa/Apj^ is the bin- integrated cross section for Higgs plus one jet as a function of pj^ divided 
by the bin width. We assume pJ^^° = 30 GeV. If we define the lower boundary of the pj, ~ 0{mH) 
region by pj. = mi{j2^ we can estimate the contribution from pj. ~ 0{mi{) to be around 30%. 

theorem, we assumed that the signal jet is of order itih- From Fig. [2| we see that this 
configuration contributes a non-neghgible fraction of the experimentally-interesting total 
cross section for ~ 30 GeV and pj. > p^*^*". Our factorization theorem holds for 
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pjKto ^ pJ^ ]^-^J^^ breaks down when pj. ~ pl^^" ^ ttih- Additional large logarithms of 

the form In^mu/PT and L x p1^^° /p'^ are not resummed in our formalism. We describe these 
terms only as well as a fixed-order NLO calculation. A different effective theory is needed 
for this regime to correctly sum the large logarithms. We do not consider this theory in 
this manuscript; our goal here is to consistently apply the currently available formalism at 
NLL' + NLO to see to what extent we can reduce the theoretical uncertainty. 

Interestingly, the p^ ~ itlh region contributes roughly 50% of the uncertainty in the one- 
jet bin, larger than might be expected from Fig. [2j We show this by computing the NLO 
cross section for an example parameter choice. We set itih = 126 GeV and p^^" = 25 GeV, 
and divide the Higgs plus one-jet cross section, whose inclusive value is ctJ^lq = 5.75;!;|gg pb, 
into two bins: the first with p;f < mH/2, and the second with p;f > niHl'^- As explained in 



detail later in Sec. IV C , we use the fixed-order cross section in the first bin since our effective- 
theory analysis does not hold, and turn on resummation in the second bin. Computing the 
cross section at NLO in each bin, and estimating the uncertainties as described in detail in 



Sec. IV C, we find 



^n1o(Pt < ^h/2) = 4.74^1:^1 pb, 

o]iM>^Hl2) = 1.0ll?|?pb. (25) 

The central values have been obtained using the scale choice for /i = nxul'^^ Although it 
accounts for less than 25% of the cross section, the region where our effective-theory analysis 
can improve the uncertainties contributes roughly half of the error in the full one-jet bin. 

We briefly comment here on non-global logarithms [48J that first occur at the NLL' level. 
Although they are not included in our current factorization theorem, to estimate their nu- 
merical effect we use the large- iVc resummation of these terms derived in Ref. |1H]- We 
include them as a multiplicative correction to our factorization formula. Their numerical 
effect is small, at or below one percent of the total exclusive Higgs plus one-jet production 
rate for the relevant values of mu and pV^^° . To check the robustness of this result we vary 
the hard scale appearing in these corrections by a factor of two around their nominal value 
of vTiu^ and find similarly small corrections. We therefore believe that it is numerically safe 
to neglect these terms in our NLL' result, although they should be further investigated in 
the future. 

^ We note that using a larger central scale choice leads to the same conclusions regarding the relative 
uncertainties of the two bins. 
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C. Scale choices and uncertainty estimation 



Since the resummation holds only for ~ mj^, we wish to turn it off and recover 
the fixed-order NLO result as pj^ becomes small. To do so, we note that the fixed order 
cross section ctnlo and the expanded NLL' cross section ct^'^l' depend only on the scale 
fJ^R = fJ'F = while ctnll' also depends on the scales fin, fJ'j, fJ^B, fJ^s, and us that 
appear in the hard, jet, beam and soft functions. The optimal choice for each scale can be 
determined by minimizing the higher order corrections to each separate component. These 
functions are then RG-evolved to the common scale fi. Consequently, the resummation can 
be turned off by setting all scales to /i , so that the full NLL' + NLO result reduces to the 
NLO one. We adopt a conservative scheme to turn off the resummation as soon as possible, 
as suggested in Ref. [SSI- In the region where ^ p^^°, we keep the resummation on. 
When p:^ ~ p^^", we switch off the resummation by setting all scales to fi, which leads to 
the fixed-order prediction. We interpolate between these two regions smoothly using 

/if*- = fi + - fi) [1 + tanh (p^ - poff) ) ] /2 , (26) 

where the index i runs over all appearing scales. We use similar expressions for the z/'s. Our 
numerical predictions are obtained using the /i*"*' expressions in our code. When pj. < pos, 
the resummation starts to vanish. We set Pos = m.ax{2pl^^° , ^) to be the default value. ^ The 
slope K controls how smoothly we turn off the resummation. We find that the interpolated 
cross section is insensitive to the choice of k. The functional forms of the interpolation in 



Eq. (26) for various values k are shown in Fig. |3j When making uncertainty estimations, 
we vary each scale separately. In the resummation region, the cross section is relatively 
insensitive to the variation of fi. In the fixed-order range, it is insensitive to fii and z/j. 

We describe here in detail how we estimate the uncertainties in both the fixed-order and 
RG-improved results. We vary all scales appearing in the cross section around their central 
values by factors of two in both directions in order to estimate the theoretical error. To avoid 
an underestimate of the uncertainty of the fixed-order calculation, we follow the procedure 
suggested by Stewart and Tackmann [32]. We split the exclusive one-jet cross section into 
the difference of one-jet inclusive and two-jet inclusive results: 

<Tij = a>ij - cr>2j . (27) 

^ The reason for this choice is that our EFT is vaUd when pj. is focated in the hard domain whose lower 
boundary is estimated to be , and it entirely breaks down when pj. falls into the "soft" regime whose 
upper boundary is roughly 2p'!^*°. 
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FIG. 3: Shown is the interpolation between the resummed result and the fixed-order one proposed 



in Eq. (26). The blue solid and dotted lines are for k = 0.04 and the red dashed and dot-dashed 
lines are for k = 0.2. When pj, is less than pos, all the scales merge to and the cross section takes 
its NLO value. We have demonstrated this behavior for the hard and soft scales in the plot. 

We estimate the scale uncertainty for each piece separately and add them in quadrature to 
obtain the scale uncertainty for the exclusive cross section: 

^li,NLO — '^>li,NLO + '^>2j,NLO • (28) 

For the NLL' + NLO result, the uncertainty is derived by adding in quadrature the separate 
variations of all scales which enter 1321 : 



^Ij — ^non-singular,/x + "^NLL'.^t + ^NLL'.^t^ + ^NLL',/xj + ^NLL',fiB,UB + ^NLL',^is,us ■ (^9) 



Before continuing we comment briefly on the structure of Eq. (29). In order to perform the 



matching to fixed order in Eq. (22), we RG-evolve the NLL' result so that all scales are set 
to the common scale fi. We then add on the non-singular NLO terms via the difference 
between the full NLO cross section and the expanded NLL' results. This explains the first 
two contributions to the above equations. As the hard and jet functions live at the scales 
^/rriHPx and p'^R respectively [H], these scale variations are treated as uncorrelated in 
Eq. (29). Finally, the variations of beam and soft functions, which live at the scale p^^°, are 



added to this. 

When we apply this formalism and assume actual LHC kinematic cuts, a large fraction 
of the cross section comes from the low-p^ regime where p^ < pos, and the fixed-order 
calculation dominates. In this situation, we split the cross section into two regions, one with 
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pveto ^ ^ and the other with > pos- For the former region, we use Eq. (28) to 



estimate the uncertainty and for the latter one, we utihze Eq. (29). We combine these two 



hnearly to estimate the scale dependence for the RG-improved cross section: 

^IjiPT > pT") = 5li,NLo(PT < Pos) + 5li(PT > Pos) ■ (30) 

Since the resummation in the result used for > pofr is turned off quickly by using the 



interpolation in Eq. (26), and the uncertainty of the fixed-order cross section used for p^ < 
Pos is obtained using the ST prescription, we believe that this leads to a very conservative 
estimate of the theoretical error after performing our RG-improvement. 

D. Numerics for the LHC 

We now present predictions and uncertainty estimates for use in LHC analyses. For 
the following numerical results, and those shown above, we use the MSTW 2008 parton 
distribution functions [51] at NLO. We assume an 8 TeV LHC, and niH = 126 GeV unless 
stated otherwise. We demand that the leading jet be produced with rapidity \yj\ < 4.5, and 
veto all other jets with p^ > p^^° over the entire rapidity range. The following central values 
are used for the scales which appear: 



H= ^J {mJj)^in{PT) min, f^H = \/ mjjp^j. , 

/ij = p^R , ^g = = Pt^° , 

'yB^^,=Xa,bV~s, ys = PT'"- (31) 



where = y m'jj- We note that these central scale values, as well as the variations 



up and down by a factor of two, are used as the /i, on the right-hand side of Eq. (26). The 
actual numerical scale choices used in the code are the /x*"*' appearing on the left-hand side 



of Eq. (26). We use n = 0.2 to produce all numerical results, although we have checked that 
their dependence on k is negligible. 

We begin by showing in Fig. [ijthe NLO and NLL' + NLO results for p'^ > 120 GeV as a 
function of p^^°, to make clear the improvement gained by adding the resummation in the 
Pj, ~ rriH region. After resummation, the scale dependence has been dramatically reduced 
for small pJ^^". For large p'!^^° the RG-improved cross section tends toward the fixed-order 
result, as desired. The pathological behavior of the fixed-order cross section for low p^^° is 
clear, as the central value becomes negative for p^^° ^15 GeV. This pathology is removed 
in the NLL' + NLO result. 

We continue by showing in Fig. the cross section as a function of the lower cut on 



for p^ = 30 GeV. Even for values of the lower pj^ cut near , a sizeable reduction of the 
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FIG. 4: Presented here are the NLO v.s. NLL' + NLO integrated cross sections with pj, > 120 GeV. 
The blue solid hne is for the RG-improved cross section and the yellow dashed line is the NLO 



prediction . The narrow blue band is obtained using Eq. (29) for the uncertainty after resummation 



while the wide yellow band comes from using Eq. (28) for the fixed-order uncertainty. 



uncertainty occurs when the NLL' + NLO result is used. The reason for this is discussed in 
Sec. IV B roughly half of the uncertainty comes from the high-pj. region, which is exactly 



the parameter space improved by our effective-theory description. 




FIG. 5: Shown are the NLL' + NLO (blue band) and NLO (yellow band) cross sections for fixed 
pj^^" = 30 GeV as a function of the lower cut on pj.. 

Finally, we present in Table |T] numerical results for both the cross sections and the fraction 
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mn (GeV) 


pJ^^*° (GeV) 


o-NLO (pb) 


'^NLL'+NLO (pb) 


A3 
J NLO 


J 1S.TT T ' 1 MT 


124 


25 


5.92+35% 

"-'•'^^-46% 


5.62+29% 

-30% 


0.299+3^% 

-49% 


0.283+33% 

-34% 


125 


25 


cr oc;+34% 
-46% 


cr c;c;+29% 
-30% 


0.300+3^^° 

-49% 


0.284+33^° 

-33% 


126 


25 


cr 7cr+35% 


cr 47+30% 
"^•^'-30% 


3on+38% 


0.284+34% 


124 


30 


c; OCT +31% 


A Oq + 29% 


n 9fic;+35% 
"•^0^-43% 


n 244+33% 


125 


30 


5 19+32% 

^•J^y-41% 


A 77+30% 
-29% 


266+35% 


n 244+33% 
^•^^^-33% 


126 


30 


cr 9+32% 


A 79+30% 
^•'^-29% 


966+35% 


n 94fi+33% 



TABLE I: Shown are the central values and uncertainties for the NLO cross section, the resummed 
cross section, and the event fractions in the one-jet bin using both the fixed-order and the resummed 
results. Numbers are given for several Higgs masses and for pV^^° = 25, 30 GeV. 



of events in the one-jet bin, f^K We define the event fraction as 



fl^ = — , (32) 



cr, 



0",;, 



where x denotes either the NLO or the NLL'-I-NLO cross section in the one-jet bin. We note 
that our values for f]^io consistent with those obtained by the ATLAS collaboration [SI] , 
which provides a cross-check of our results. The total cross section amc, as well as its estimated 
uncertainty, is taken from the LHC Higgs cross section working group [55j . The uncertainties 



shown are calculated as discussed in Sec. |IV C[ Results are given for mu = 124 — 126 GeV, 
and for p1^^° = 25 and 30 GeV. The reductions of the uncertainties are significant for both 
values of pl^*^". Symmetrizing the error for this discussion, the estimated uncertainty on the 
cross section improves from ±40% at NLO to ±30% at NLL' ± NLO, a reduction of one 
quarter of the initial value. The one-jet fraction uncertainty decreases from ±44% to ±34%. 
For pJ^^° = 30 GeV, the error on the cross section decreases from ±36% to ±29% when 
resummation is included, while the error on f^^ decreases from ±39% to ±33%. Numerical 
results for other parameter choices are available from the authors upon request. We note 



that these are extremely conservative error estimates, as discussed in Sec. |IVC[ We default 
to the ST prescription over a large region of the relevant parameter space, and turn off the 
resummation at a relatively high value of pj.. Enough of the error comes from the high pj, 
region that our RG-improvement is effective in taming the uncertainty. 

V. CONCLUSIONS 

We have studied in detail the resummation of a class of large Sudakov logarithms appear- 
ing in the perturbative expansion for Higgs production in the one-jet bin. These logarithms 
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occur when the jet transverse momentum is of order the Higgs mass. For certain param- 
eter choices they lead to pathological behavior of the fixed-order perturbative expansions, 
including negative cross sections and scale-uncertainty estimates that do no properly ac- 
count for missing higher-order corrections. Past attempts to handle this problem increased 
the theoretical error estimate. While theoretically correct, this has the unfortunate effect 
of introducing a large systematic error into experimental analyses. The results we present 
here tame the poor behavior of fixed-order QCD by controlling the jet-veto logarithms to 
all orders in perturbation theory, leading to a more reliable theoretical prediction and a 
reduced uncertainty estimate. We have reviewed the necessary formalism to understand the 
resummation, and have extended the theoretical accuracy to the NLL' level by calculating 
the full one-loop soft function for this process. Using our matched NLL' + NLO result for 
the cross section, we have performed a detailed numerical study of exclusive Higgs plus 
one-jet production at the LHC, including an estimation of the non-global logarithmic effect 
and a careful accounting of the theoretical uncertainties before and after resummation. The 
estimated theoretical uncertainties in the one-jet bin are reduced by up to 25% using the 
NLL' -|- NLO prediction. This is an extremely conservative error estimate, as argued in the 
main text. 

We believe that it is now time to revisit the theoretical error treatment used by the LHC 
experiments in their studies of Higgs properties. Current studies use fixed-order perturbation 
theory, and the error treatment suggested in Ref. [32], to estimate the uncertainties induced 
by dividing the signal into exclusive jet multiplicites. While this is the best choice to correctly 
handle the uncertainties induced by the jet veto if only fixed-order predictions are available, 
the resummation of jet-veto logarithms is now known for both the zero-jet and one-jet bins. 
These predictions do not exhibit the pathological scale variation that led to the prescription 
of Ref. [32] • We are confident that this is only the first of many upcoming improvements 
for predictions of Higgs-boson production in association with a fixed number of jets. A new 
factorization theorem for the low-p^ region for exclusive Higgs plus one-jet production, an 
improved treatment of In R effects, and new high-precision fixed-order calculations will all 
further reduce the theoretical uncertainties and help provide a sharpened image of the newly- 
discovered scalar at the LHC. We encourage the experimental communities to incorporate 
this new knowledge into their analyses. 
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Appendix A: Fixed-order jet and beam functions 



In this Appendix, we tabulate all ingredients needed for resummation at NLL' accuracy. 
We start with the NLO calculation of the jet and the beam functions, whose operator defini- 
tions can be found in Ref. [H]. The anti-^T jet function is calculated using the measurement 
function 



M., 



e(Ar/^,+A0?.<i?^) + O(pr 



(Al) 



We explicitly check that the coUinear radiation leaking outside the jet is power-suppressed 
by p^*° after correctly subtracting the soft zero-bin contributions. Since numerically R <^ 1, 
we can simplify the measure using 



Ar^J + A02. = 2cosh(A77,j) - 2cos(A0y) + 0{R^). 
In this limit, the NLO jet functions for gluons and quarks become 



(A2) 



-^9 



27r 
2tt 



67 _ Stt^ 
13 Stt^ 



23 Uf ^ , u 

^ +/3o\og-f 



4 



+ 3 log 



9 2 



+ 21og^ ^ 



+ 0{R^) 



(A3) 



p^R p^R 

similar to what was obtained in Ref. [52] using a slightly different jet algorithm. We note 
that the jet functions are normalized so that the leading-order results are unity for both 
quarks and gluons: jf^ = 1. 

The measurement operator for the beam function with a single emission is 

MB = e < pI''") e {\r]i\ < r^cut) + e{h\ > Vcm) . (A4) 

Experimentally, r/cut ~ 4.5. For simplicity, we set rjcnt = oo here. We note that this difference 
does not affect the anomalous dimension of the beam function, it only changes the finite 
part. The calculation is performed using the 't Hooft-Veltmann scheme. The NLO matching 
coefficient X for the various beam functions are found to be 



27T 

27r 
2t\ 
27r 



Pt 



4 log log ^5(1 

Pt 



-2pgg{z) log 



n-p 

FT 



-2pgg{z) lOg-^ 

Pt 



2pgg{z) log 

■ 2^99(2) log 



Pt 



fx 



FT 



+ Z 



2z{l-z] 



(A5) 



20 



with 



2z , , \ — z 



, , 1 + fl -2 



z 

2 



Pgg{z) = 1 - 2z + 2z\ (A6) 

We note that the beam functions are normahzed so that the diagonal matching coefficients 
are delta-functions at tree-level, while the off-diagonal ones vanish: X-^°^ = 6{1 — z), xf^^ = 
for i ^ j. 



Appendix B: Anomalous dimensions and RG evolution 

The beam, jet, soft and hard functions appearing in the factorization theorem of Eq. ^ 
all satisfy evolution equations of the form 

= r^(/^)^(/^) • (Bi) 

The soft and beam functions, which also contain rapidity divergences, have similar evolution 
equations in the rapidity regulator: 

These are easily extracted form the poles of the one-loop calculations for each of these 
objects. The general solution to these RG equations can be formally written as 

u) = U (yjj,, u, /io, t^o)F{hq, uo) . (B3) 

These RG-improved expressions for the beam, soft, jet and hard functions are then used in 
Eq. The initial conditions in Eq. (B3 ) can be determined from the fixed-order calculation 



presented in the previous section and in the main text for the soft function. 

From these fixed order calculations, we can determine the anomalous dimensions used for 
RG evolution. The anomalous dimension for the jet function is given by 

r,^ = 2r,,,p7;2 log ^ + tj, • (B4) 
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Here, i labels the parton flavor, and can take on the values i = q,g. For the beam function, 
the anomalous dimensions are 



Tb = 2T^us,T^ log , 
Pt 

= 2r,,3p7;'iog — + (bs) 

n-p 

where = Ca for gluon and Tf = Cp for quark. The anomalous dimensions can be 
obtained as expansions in the strong coupling constant. For the resummation performed 
here, we need the following terms in each expansion: 

as /as\^-p 

Tcusp = ^ro+(^-j ri + ..., 

'-^s i '^s i fT-)r\ 

IB. = ^7o, 7J. = ^7o- (B6) 

We note that the non-cusp anomalous dimensions of the beam and jet functions are the 
same. We have the following expressions for the necessary anomalous dimensions, as well as 
the relevant coefficients of the QCD beta functions needed both here and later: 

To = 4, 
ri = 4 

7o' = 2/3o, 7o' = 6C^. (B7) 
The RG evolution of the soft function can be easily obtained by direct differentiation of 



Eq. (20). As the soft function is sensitive to both the jet and beam directions, its anomalous 



dimensions take on a more complicated form than those of the other quantities. We find 
= 2r_p { {T! + T^) In^ + yj{T^ -Tj-T,- Tj) ln-|^ | , 

= -2r™sp (T,^ + T,2) In-^, (B8) 

Pt 

where a, h denote the beam directions and J the flnal-state jet direction. We note that 
+ F^^ + F^^ = 0, as required by RG-invariance of the cross section under //-variations. 
The anomalous dimension of the hard function was studied in detail in Ref. E3], and we do 



not reproduce it here. It can be derived from the above expression using RG-invariance of 
the cross section: F|^ + F^^ + F^^ + F j + F^ = 0. 
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Appendix C: Solutions of the RG equations 



We reproduce here the solutions to the RG equations presented in the previous section, 
which are required in Eq. (B3). The evolutions of the jet and the beam functions are given 
by 

/ \ -2T]Ar(tMJ,^l) 

UJ,{HJ,^J^) = exp[-2T^S{fij,fi) - AjXfij,fi)] 



UB,a(,f^B,l^B,fJ^,l^) = exp 



,2 1 



exp 



-^a ^r(/^s, /^) log ^ - (/xs, /x) 



(CI) 



The solution to the RG equation for the hard function is is 
UH{fJ'H,fJ') = exp 

i 

2TiAr{fiH,tJ.) 



X 



n 



(C2) 



where we have set ARj^ = e Ai?jj = e'^-' and Ai?^^ = 1. We also set Ui = if i = J, 
otherwise we have Ua = Xa^/s. The soft-function evolution factor is 



Us{fis,i^s,fJ^,i^) = exp 



T,; ■ T, 



(C3) 



For the NLL resummation, we need the following factors: 



Ar{fii,fif) 



2/3o 



logr + 



47r Iro /3o 



fr- 11 



(C4) 



and 



47r 



- - logr ) + ( ^ - ^ ) (1 -r + logr) + ^^logVl 



(C5) 



where r = as{\if)las{\ii]. Aj/b, Ah and are needed at leading order, and can be obtained 
by substituting the Fq in Ap with the corresponding 7q and expanding in a^. 
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